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A DESIGN PROGRAH FOR 


SUPERCONDUCTING ELECTRICAL MACHINES 


by 
DAVID GREENEISEN 


ABSTRACT 


Submitted to the Department of Naval Architecture 
and Marine Engineering on May 17, 1968 in partial ful- 
fillment of the requirements for the Master of Science 
Degree in Electrical Engineering and the Professional 
Degree, Naval Engineer, 


This paper presents a procedure for the design of 
superconducting electrical machines. The magneto- 
static field problem of a cylindrical superconducting 
machine is solved. Equations for the magnetic fields 
are developed, and from these, oxprossions for the 
various machine parameters are obtained. These expres- 
sions are adapted to a computer solution of the design 

| problem. The computer program obtains the design 
parameters for & minimum volume machine dosign. 


Peces1en Study for a typical marine clechriucal 
propulsion system is conducted,- The results of Chas 
study indicate quantitatively the reduction in weight 
and space to be obtained in using supcorconäucting 
electrical machines vice conventional electrical 
machinery. The study aiso presents some general 
characteristics of superconducting machines, 


Thesis Supervisor: Н. Н. Woodson 
рее Professor of Eloctrical Enzinsering 
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I. INTRODUCTION 


A. Background 
The relative flexibility and simplicity accruing to 


electrical power systems, in comparison to mechanical power 
systems, has long been recognized. Specifically, electrical 
power, and transmission systems provide flexibility of 
installation, are easily suited to automation and have a high 
degree of dependability. 

Dospite these advantages, and others that could be obtained, 
the full benefit of electrical power systems has not been realized 
in marine propulsion applications. To be sure, electric drive 
systems have been used in ship propulsion, however, such use has 
been mostly limited to small or intermediate size propulsion 
plants. The primary reason that electric drive systems have not 
been adopted for large power plants has been that the weight and 
space for such systems is greater than that required for geared 
propulsion systems. In addition, increased cost (as much as 20%) 
and lower efficiency (85% for DC) have been contributory Da), 

It is advantageous to consider some of the benefits that 
could be obtainsd with electric propulsion in order that aspects 
eats UE lization micht bo jJuaged in proper perspective, 

As mentioned above, an electric propulsion system affords 


As 


greater flexibility of installation. In constrast with a 
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conventional goared propulsicn system, extensive shafting is not 
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required and the location of the power source is not as 
mestricted. In addition, the large, noisy roduction gear could 
be eliminated. The noise generated by this component is cf 
increasing military significance. 

If a cycloconvertor^ is used, additional benefats derive. 
It would be possible to use a constant frequency generator, 
thus permitting a more efficient design of the prime-mover. 
This will also permit greater use of newly developed power 
sources such as the gas turbine. Added to this would be more 
adaptation of remote and automatic controls, plus much improved 
maneuvering E CC 

Thus it seems that if the disadvantages (weicht, space, and 
possibly efficiency) of electric propulsion are reduced, it 
will be given greater consideration as an alternative to geared 


propuision in larger power systems. 


B. The Role of Superconductors 
It has been realized for some tins that when the temperature 


of an electrical conductor is reduced to a few degrees absolute, 


the resistance of the conductor decroases to an unmeasurable 


—+ 


amount. The implications of this fact have considerable value 


for use in electrical machinss. If ths resistance in the windings 


of an electrical generator or motor is removed, much greater 


5 e ° ç ` Р е - > en a E S ° - 
Current densities can oe achieved with an acconpanying incvease 
= ch E 4$ 3$ A] ` m ` < eevee a T STE ° ^4 E ~ ( 3) LY LLY, 4 
in ө macnevic flvx densivies within the machines. jas opo 
sur 
a’ 

^ us \ £ بم‎ © u /7 y 2. ry 9 ya 16 Е en e „irn Bi ” = f Es A q de NN a A ^. e Ya e A Ph > : 
Si ос E o E Dome oui ar alerts mu, ao ОСИ 1 
q 


Е 


е ТЕА 4 c Pm E T OA : “ a ыа ; 
voltage and frequency. Sce гоГогөпсо ц. 





oe 


increase in flux density is obtained, then greatsr power can 

be generated within a given volume. Thus it follows that the 
volume of an electric machine can be reduced if superconducting 
windings are used. Studies to date, 13) have demonstrated the 
significant reduction in volume and weight in superconducting 
machines. 

However, the application of superconductivity is not without 
its complications. The first question which must be asked 
concerns the method of operating a machine wherein all, or at 
least some, of its parts are at a temperature in the region of 
liquid helium. 

It is at this point convonient to limit the scope of the 
examination. The ensuing discussion will treat only machines 
with a rotating field winding, operating in the superconducting 
region, enclosed by a stationary armature at ambient tempsratures. 
This particular configuration is only one of several that cculd 
be considered. 

ihere are two particular reasons uhr this configuration is 
selected. The more obvious is that it would be impractical te 
transfer the high machine power output through siip rings, eS) vould 
be necessary if the armature were rotating. The other reason 
vesults from the fact that large heat losses are encountersd 
when superconductors are subjected to slternating currents. ‘3? 


Thus, only the field winding is made superconducting. 
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Nonetnoless, most of the techniques in this paver ar: 
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With the discussion constrained to considering a rotating 
superconducting field winding, it is possible to dilineate some 
design procedure. Thus far, it has been found that the use of 


a thermal bottle, or dewar”, filled with liquid helium and 


encompassing the rotor, is a satisfactory solution to the 


problem of maintaining the vindings in & superconducting state. 
Several techniques are under development fcr introducing the 
helium into the rotating dewar and passing electrical leads 

into the low-temperature region. Generelly it is favorable 

to use the liquid helium dewar in conjuncticn with a liauid 
nitrogen dewar inmediately exterior. The helium maintains thse 
rotor at about haar while the nitrogen offers a 770K buffer 

for improved cryogenic performance, At first, one might think 
that the refrigeration load for maintaining such lou temperatures 
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ple 


would bo excossive. However, 


s very little heat 


Ho 


the low temperaturas are obtained. thers 
generated by the windings. Indeed, heat introduced into tho 
low-temperature region through the shaft and support structure 
is more significant. А reasonable estimate of the power 
requirements for refrigeration during operation of a super- 

(1) 


conducting machine would be about 0,5 percent of rated pover. 


in addition to the requircment that a superconductor be 


03 


maintained at a low temperature, thore is a limitation on the 
amount of current carried while exposed to a magnetic field. 


As indicated in Figure 2, if the megnetíc flux density increases, 
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provided to confine the magnetic fields within the machine. 


ас 


Ehe регті зи 01е current density is reduced. If current or 
field exceed the limits imposed (the so-called critical 
values), tho conductor reverts to its normal state. Thus 
Anis necessary, in any particular machines design, to 
determine that combination of maximum flux density and 


current density which best suit the design. 


critical current 
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[шо Corrente rux Relationship: Iorio 
Typical Superconductor 


The very fact that a superconducting machine develops 
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extremely intense magnetic fielós invitos yoei 
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complication. Since fields of the order of several kilo- 
gauss can be expected within the machine, disturbance of 
equipment in the vicinity of the machine wiil certainly 
result, as well as unbalanced loads on the superconducting 


field winding. Therefore, it is necessery that a shield be 


( 
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This shield would typically be of laminated steel con- 
struction and will be a Significant portion of the size 


and weight of the machine. 


oU. Design Proeram 





Though superconductins machines of the size susgested 
by the foregoing discussion have not yet been built, they 
are certainly feasible. Anticipating the construction of 
such machines, it is considered desirable to investigste 
some of their expected properties. Since there is vested 
interest in the reduced size and weight of these machines, 
it is natural that a desisn procedure be developed for 
obtaining minimum weight and volume. 

ane objective herein is to find a procedure for finding 
design characteristics of least size machines, The procedure 
developed is founded on a field-theory solution of the 


electric machine and is adapted to a computer solution. 


272 
II. PROCEDURE 


A. Development of a Model 

The first step in the analysis was to obtain a 
mathomatical model for the magnetic fields involved. This 
model was necessary to determine the parameters which would 
be used in machine design. As mentioned previously, the model 
and attendant equations are for a rotating field winding inside 
a three-phase stator winding. 

Analysis proceeded by considering both the field and stator 
windings to consist of cylindrical current sheets. Equations 
for the magnetic fields generated by each winding were develcpoa 
separately and later combined to obtain necessary inductance end 
energy expressions. 

To illustrate the analysis, consider the field windings. 

As illustrated in Figure 3, a cylindrical current sheet, K(r). 


of thickness dp , located at somo radius r, R. ¿rno is examined 
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Fourier analysis of this current sheet yields: 


J ar 





sinnpg. (II-1) 


Knowing that these current sheets exist between r * Ry and 


p =R and that the fields associated with each sheet ars 


2? 
Laplacian, enables determination of the field expressions. 


ver finding an expression for the magnetic field for r : Но, 


(5 


ons can find the contributicn due to the shield. Since at 


r Een. with an infinitely permeable shield, 


ne 
SO 6 
oi, au © = O, 


H = finite (11-3) 


an expression is found for the fields induced by the shield. 

The expression thus obtained is then added to that for the 

Fields generated by tro field windin 
nis same procodurr obteins Tor detemminine the magnetic 
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fieid system resultin from the armaturo windin 








Figure li Distribution of Armature Current Sheet 
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атшайпгө Windings are prosentod in Tables l and 2. 
Using tho equations for the magnetic fields, an expression 


for encrgy stored in tne armature field system is obtained from 


et 5 
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Equating the result in 


N E Ta 
Wa 7 > L la: (11-5) 
where 
2 ‚2 
pow, RI Р (11-7) 
6 H. 


yields an expression for the armature self-inductance. A 
similar calculation is performed to obtain tho field winding 
self-inductance. 


Mutual inductance is found from 
DEN =) , (11-8) 


where A is the total field-winding flux linkage of the armature 
e. a a s 2 

windings. By recosnizing that ӨМ Ит is - К.) expresses tho 

number of armature turns per unit cross-section, then the 


second spatial differential of the flux linkage is expressed 


by 
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Imtegsrating this expression over the cross-section of the 
armature yields the total flux linkage and thus the mutual 
inductance., Expressions for field and armature self- 


inductance and mutual inductance are presented in Table 3. 


B. Design Equations 


If one considers the amplitudes of the line-to-neutral 


armature voltage and armature current, one obtains 


V. “WT, (11-10) 
2 2 
І = т (0 ii R;)J, (TI) 
a EE Bie 9 
a 
DESI (11-12) 
£ oan a” 5 
where 
2 2 
I, = T(R5 - Ry)Je (11-13) 


2N a ° 
í 
This power is internally gonorated power, which must 


be reduced, to account for armature impedance and tho load 


power factor, in order to obtain rated output power. 





3$ s ME id [Ip ci T ie * Pe dc^ d mS 
lucc a Vector Diagram of Generator Volteges 
m E nee pres. m 





-16- 


TABLE 3 
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From the vector diagram in Figure 5, where 


La = Dita s (11-11) 
the ratio V/V, is obtained. Then rated power is 
P= ор cos ё. ‚ (II-15) 


In developing the equations describing the machine, sone 


approximation was necessary to account for the effects of the 


(10) 


end-turns. Figure 6 shows the geometry involved in making 


this approximation. | 





Figure 6 End-Turn Geometry 
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of this coil is Im/3p radians. Since this coil is at a radius 
r= d(R + R,) (11-16) 
со i 
the approximation can be made that 


aisi л ER + Ry) = (Ry + R,). ОТ 


(6) 


Using this expression, it has been found that a reasonable 
length to assume for self-inductance calculations is the total 
length 

DM Mr ONIS, (11-18) 


whereas the length to be used for mutual inductance is 


а + dee (11 195) 


C. Design Limitations 

As described in the Introduction, the low resistance 
Character Of a Superconductine coil is not only a TLunstierzor 
temperature, but aiso depends on its current density and the 
Magnetic field to which it is exposed. Thus for a given design 
it is necessary to determine the maxi:nr magnetic flux density 
associated with any particular current density, to ensure that 
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limitation are detailed in the program construction techniaues 
ei Appendix C, 

As also mentloned in the Introduction, the machine design 
considered employs a shield (constructed of laminated steel) 
which contains the dense magnetic fields within the machine. 
For this shield to be effective, the fields it constrains must 


be less than those which will cause saturation of the shield. 


If the shield material saturates at some density B rated! then 

ono limitation is 
B _ (=R) 5 (11-20) 
max. | O rated” 


Another limitation exists in the requirement that the maximum 


flux density carried by the shield must never exceed B 


rated 


As illustrated in 


I 


PXeupro 7/  Shiieia Flux -Patcern 
Figure 7, tho total flux accumulated by the shield over an ansle 


of "/2p must not be sufficient to cause saturation in the cross- 
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section at A-A. Since 


w/op 


Bei, = ; Bg (Ry'1R,d9 , (II-21) 


a final constraint is obtained in 


Z . 
ф (R - СКК (ТТ=22) 


nu LS 


There are two final considerations which may govern a 
machine design. First, the heat dissipation in the armature 
winding must be considered. Generally the maximum permissible 
current will be determined from the linitations imposed by the 
method adopted for cooling the armaturs. For purposes of this 
design it was assumed that an average armature current density 
br 1.6 x 109 amp/m* could be accepted with chilled water cooling 
of the winding. Finally, one must consider the speed of the 
machine rotating parts. The strength of the rotor structure 
mill provide a limit to the rotor tip speed that can be sustained, 
On consideration of previous conventional machine design it was 
determined that 700 ft/sec was a reasonable upper value for tip 
Speed. However, since methods of constructing superconducting 
MAndanes are far from beine fixed, this value гог шектш р 
Seeca 1S spoculatave. The tip speed of the rotor will also 36 
governed by the performance of the rotating dewar. In this 
respect, the Pluidic and thermodynamic properties of liauld 


ШО л, under hiosh Centriiuge: forcec, wl) goyeen. neuste 
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aspects of this considoration are yet to be determinea, it is 
assumed here that the previously stated maximum tip speed is 


acceptable. 
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III. RESULTS 


The significance of this work lies in the results of a 
design study which was conducted, employing the design 
procedure described in Part II. The design study involved 
what was considered a typical marine electrical propulsion 
system; a turbino driven generator rated at 26 M7 (35,000 hp) 
at 3600 rpm, and two direct drive propulsion motors rated at 


17,500 bhp (13 MW) at 300 rpm. Other paramnoters assumed wore: 


Je = Je NX 109 amp/m“ 
Jo = еб. K 109 emp/m* 
= 5.0 web/n^ 
peas 20 web/m 
о, = 1000 ke/m> 
0, = 1000 kg/m? 
(copper with .li5 packing factor) 
р. = 8000 kg/m? 
(steel) 
§ = . OEM 


Tip spoeds of up to 21h m/sec (700 ft/sce) wers considered 
in both the generator and motor designs. Also, for the 
generator, designs with power factors of 1.0, .850 and .707 


wore examined. Desirn characteristics of machines wi 
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The computer results and graphs of the characteristics 
of some of the more interesting designs are included in 
the following pages. 

One quantity, not calculated by the computer progran, 
thet should be examined is the machine efficiency. 

Since 


к Erat ТТЕ, 


х 
where ps Tosistivity м От 10-8 ohm-m and cross-section 


area is 
2 2 e і 
Nu ui (III-2) 


A 


for a .uS packing factor, then for the optimum generator 


Besten (pf = .707, p = 6, = 21l, m/sec) 


R,/NÉ = 6.7 x 1079 oh turns. 
Enowins N Ia?’ 20-195 рост лелро гла 


IR. = 965 watts. 
a 8 
These losses, however, are for only one pnase and were 
computed with a peak vaiue for I, anstead оГ тиз, Theor CIO 


resistive lcsses for all three phases are 


ту 965 * E Seo aos e 


af an equal amount is allowed for eddy current losses, a 
figure of 2990 watts is obtained for the total losses. 
since this ‘Is for a 26 lil! generator, this loss will cause 
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Figure 8 
Generator Specific Weight 
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Figure 9 
Motor Specific Weight 
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IV. DISCUSSION OF RESULTS 


ile moso significant result cí tnis study lies an the 
low values of specific weignt (lb/hp) for superconducting 
machines. As indicated in Figure 8 specific woights on the 


m^ 
es 


order of 0.25 lb/h> ` are theoretically possible in machines 
of the size typically founá in merine propulsion plants. 

The results of the study of direct-drive propulsion 
motors were less conclusive. It is apparent from Figure 9 
that specific weights of 1.40 lb/hp are possible and that 
much lower specific weights are likely to be achieved in 
machines with larger numbers of poles and higher tip speeds. 
Nonetheless a specific woight of 1.4.0 is still an improvement 
over a conventional machine. This is evident from a 


comparison with the data for the propulsion motcr in the 


USS Hunloy, AS-31, given in Table L. 


hp 15,000 
pole pairs 26 
rpm Dom 
1b/np 10.9 


ç 


Table Ë Characteristics of Propulsion Kotor 
(11 


T 
OL USS HUNI st 


Consideration of the results of the generator designs 
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bo a limit in machine design. However, this is not likely 
to be the case in lower speed machines such as the direct 
drive motcr. 

Machine characteristics such as low X, and high efficiency 
are expected as a result of the smaller amount of armature 
winding required. 

Aside from the significance of the data given for the 
various machine designs, there were two other notable 
Besults of this work. As described in the procedure for 
determining the field winding design (detailed in Appendix B), 
in order to reduce the number of variables involved it was 
assumed that the fields induced by the shield did not affect 


B However, the study conducted in Appendix I indicates 


pmax” 
that in some instances (especially for small values of p; 
this assumption is not warranted. From Figure 16 it is 
found that for the optimum generator design (x = .822, 
y = .848, p = 6) the error introduced is about one percent, 
whereas for the optimum moto» (x 7.801, y = -&87, p = 6) 
the error is about .80 percent. These accuracios aro ox- 
pected to be much greator than the accuracy of the math- 
ematical modsl being used. 

Another result that became evident during the design 
study was that the linear interpolation techniocue used to 
find the minimum volume design (described in Appendix D) 


was nob sufisclenely accurato 312 all 3ngtences. Since cle 


MINIMA OL se volume vs. R. би уо je GULo pron Lo» 
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D —1l and 2, the linear interpolation yielded erroneous 
results in a few designs. 

The accuracy of solution provided by the design program 
(within the limits imposed by the assumptions as discussod 
ín Appendicos B and E) is arbitrary. In the design study, 
a variation of 1% in output power was allowed. Likewise, 
in instances of solving for maximum flux density, the rate 
of change was required to be less than 107 at the calculated 
maximum. By changing these tolerences, the designer can 


obtain the accuracy desired in the solution. 





SC 
V. CONCLUSIONS AND RECOMMENDATIONS 


From the results of this study, the following can be 
concluded: 

1. Superconducting electrical machines should provide 
considerable savings of weight and space in electrical 
propulsion systems. 

2. The electrical efficiency of superconducting machines 
is substantially higher than the efficiency of conventional 
machines. 

3. The solution of the design problem for a super- 
conducting machine can be mace at least as accurste as the 
design of a conventional machine. 

l. High speed machines will be limited in speed of 
rotation by the maximun tip speed that can be sustained by 
the rotating parts. 

5. Ап optimum design for a slow speed machine will 
roquiro a largo number of poles, 

in the conduct of further investigation of the optimal 
design problem, or in application of the design program 
presented here, the following is recommended: 

l. The body of data (array Blp, y)) representing the 


graph of B const, vs, y, Should be expanded to account 


P mnax 
for larger values of p. 

с» The design program should be modified to allow 
variation in VIIP. Аз indicated by Figures 8 and 9, an 


JAN 


ОР ОЛЕ л лусу Lor somo valuo ol WEN 
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that is less than the maximum VIIP, By searching over 
a range of this variabie, & numerical techniguo could 
easily be used to find the optimum design. 

б, When a design is obtained, 1t is necessary Lo 


ensure that the error in B as described in 


pP max’ 
Figures 13 to 16, is not excessive. It may be desirable 
to incorporate these figures in tho design program, as 
was done with array B(p, y). It would then be possible 
to check each design within the program, making appropriate 
adjustment in the design parameters until a satisfactory 
result is obtained. 

H. To ensure greater accuracy in the subroutine MIN, 
described in Appendix D, a false position solution ne 
should be substituted for the linear interpolation method. 


A basic routine, such as the one presented on the following 


page, might be considered. 
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APPENDIX A 


Input Data 

In organizing the machine program thal was ultimately 
used, it was necessary to anticipate what information was 
likely to be available or fixed for the machine design. It 
was anticipated that certainly the desired power output 
would be known. The angular velocity of the rotor and 
power factor of the load may also be determined fron 
particulars of the intended machine application. Further, 
t was expected that materials selection would precede this 
portion of a design, thzrefore various matcrial character- 
istics will be available. Also expected is some knowledge 
of construction technlauss to be employed, e.g. the type of 
armature cooling used and characteristics of the rotor 
design, Finally, some information concerning the dewar is 
necessary. Specifically the anticipated thickness of the 
dewar wall must bə provided. Thus, input data will include: 

as power output, PS 

b. mechanical frequency, Un 

С. Maximum armature current density, Je 
d. dowar wall thickness, A 
e, maximum tin speod of dewar, 


(fixed by strength 


К 2 $ 7 е 
t 3j cu. in roor doci neon 
limitations impos y 


Con 
ed by dowar te chnolcgy) 


Bone t c OL 
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g. flux density to cause shield saturation, B ated 


h. average density of rotor, armature and shield, Oo» Par Lg 


As mentioned in the Introduction, the superconducting 
material used in the field windings will have characteristic 
walues of critical current and critical flux density, Thero- 
fore, additional inout data will be a value of current density 


in the field windings, J and the associated value of critical 


T3 
miux density, В ° 
max 

For any field winding design, the geometry of the magnetic 

telas can be found. If Je is known, the maximum flux density 


can be calculated. This maximum flux density must be found 


to ensure that it does not exceed the critical flux density 


€ 


IT, E i ` 


associated with J p (B Ju 

Conversely, if de ena Bey ate given, cortan 
limitations on field winding design are imposed. 

It was found expedicnt te provides a body of data which 
could be used in the comvuter program for finding the best 
meta winding design for sone given J. and о The 
method for obtaining and utilizing this data is described 


in Appendix B. 
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APPENDIX B 


Input Data for Field Winding Design 


In the design of the field windings, three assumptions 
are made. First it is assumed that the fundamental component 
of the magnetic fields is the significant part, This assump- 
tion is fairly standard and has been shown to be acceptable. (©) 
Second, it is assumed that the fields generated by the arma- 
ture may be neglected. Since the armature fields oppose the 
fields generated by the field windings, this approximation 
is conservative., Third, it is assumed that the fields 
reflected by the shield may be neglected, Although these 
Fields add to those from the field winding, they arə smali, 
becoming insignificant as p increases, These approximations 
ere made in order to decrease the number of variables in the 
calculations. 

From Table I, it is seen that in the region Ry р Ro 
the maximum field intensity will be in the p-direction 


for 8 = 0. In accordance with the sbove approximations, 





for p F 2, 
Эл RA -p*2 H | 
Rr a „2р + (2 + р)| = SS pue p 
: 2) ei Le 
(4 -p”) | M" 
or 
| 0 +2 
| ee 
const. " >, J h (limp с Wee 
For 
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ІГ, for a given value of y, 


d(Bp /const.) И 
G 0 (B-3) 


is solved for (E/R,), and the value obtained used to 
evaluate Eq. (8-2), tho maximum value of Bg /const. is 
obtained for that y. 


For the case p = 2 





2J gl R 
f 1 1 2 
H =z —- l- me B- 
т T o п C e (B-l) 
and 
B B y TM R 3 R 
ааа = сЕ = 1. 0 у =" + e in 58 
const. 2ш ello i it, О R. 
(B-5) 
Thus, if Ro милоти, ас зироЕ Ерте Lo obtain 
E max 
Емал. z: түру (B-6) 
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In the following computer program, for the IBM 1130, 
B maz) const, is computed for values of p = 2 to 6 and 
y = 0,05 to 0.95, Tho results arc disvlayed sraphically 
WE таро О Using these resolus, 16 18 possible uo 210 
the value of y corresponding to a given B 6 сом CNG ip 


thus fixing the radial cimensions of the field windings. 


Be Coe. Ds sols NOt Somnusedzesinecgtor Ds al usa rp 
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In the program, the following notation is used: 


p = number of pole pairs 
y š Ry/Ro> Ry is tho inside field radius, etc. 
R = @/R,, i.e. Rl in the program is the first value 


of р/н, and is not the inside field radius. 


p = Bp /const., 


) 


B MAX = B o may) const, 


+ * n 


The program procedure is simply a solution of Eq. (B-3) 
and evaluation of Eq. (B-2) or Eq. (B-5) for the maximum 
value. For the case p = 2 a false position nuc is 


used to solve Eq. (8-23), wheroas the Newton-Hephson 


б) 


technique! is employed for all cases p ў 2. 


The printed results of this program include values of 


p, y end ans RE, wnich are used in the design program. 


fhe data represents the curves in Figure 10. Within the 


design program /const, is computed, with which it is 
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* е. е 


tenen possible to obtain fron tho curves a velue of y wnich 


will satisfy the constraint imposed by the current and 


Be Ого сашсизас Lo solve Tor the date an cho arras 
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B(p,y), and tme subroutine used to interpret В(р,у) to 
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find a value for y, are included in the following pages. 
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APPENDIX C 


Design Program 


The design computer program, written for the IDH 1130, 
is an iterative solution of the basic design equations to 
find that machine design which requires the least volume. 

The program uses the normalized variables x and y which 
represent the thickness of the armature and field windings. 


The variables R, and R. serve to fix the radial location 


2 , 
mM LNG windings. The machine length is found from consiaer- 
ation of the working volume required to obtain the necessary 
internally generated power. The reduction from internal 
power to rated power is obtained from Eq. (II-15). The 
itorative procedure followed is based on varying values 

of Ros which was found to be the most convenicnt parameter 
for this purpose. 


Since 


Vol, = T Bos і 


the variation of Ros with R. Was i1apbst Oxaminogo. By 


evaluating Eq. (II-21) and substituting in Ба. (11-22) опе 


obtains | rr pie , 
а 8 
es На EN (022) 
x TPPratoag 2 Р) 
o. 
dic R (0-3) 


35) Q 
ә, 
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Кор! 
DE P 
кле Variation of Ros wach R 


For small values of Ros the shield is in the vicinity 


of intense magnetic fields and R Кы must be large. Аз R. 


increases, the thickness of the shield decreasos but Ros 


FS lll larger than Re . The value of i which requires 


the smallest Ros 15770020 со 06 


> 


ES. р 
AN (pC, ) (C-L) 


Though this value of Ro does not necessarily yield a 


3 


fe 


minimum volume machine (since no variation of 1 


t 


considered), it is gensrally near the minimum volume value 


becauso of the significance of R іи Ед. (0-1). 
An R. obtained from Eq. (C-h) will satisfy the constrain? 


of Ба. (11-22). Hovever, it must also satisfy Eq. (11-20). 





Since 
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T p S o — “rated 

(C-5) 

it is required that 

DE 
R. >` (0,/p) : (C-6) 
since 
Ded = 

= > — n 
us opt. on ED (6-7) 


Eq. (11-20) is also satisfied. 
Turning now to the design program, the following 
procedure is carried out for successive values of p: 


From the input dava, one obtains 


= Bo 
bvaluating 
B aan’ erst. са (C--10) 
HN cue e Ko р 
and entering the curves of B , ./const. vs. y. described 
И С.Л. 


in Appendix B, a value of y is obtainsd uba Е 
information it is possible Lo evaluate ©. anc Find R 


which is used to calculate an initial mechino desi 





-52- 


At this point in the program, the difference between tho 
internally generated power and the desired power output is 
unknown. Assuming that this difference is largely due to 


load power factor, it is approximated that 


РЕ = P_/pf (Cane) 
Since 
P. us (6-13) 
= 2 Ww MILI (0-11) 
2 e 1a 
where 
Qs * pO, (C-15) 
| SN 
م1‎ = п(2-у ) 5 Je (0-16) 
2» per 
ПОЕ 2 
ON, 


and the appropriate expression for M is selectod from 


Table Ii, one obtains 
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mor p = 2, From Eq. (C-18) or Eq. (C-19), it is possible to 
solve for 1. Then, in accordance with Ea, (11-18) aná 


ШО. (11-19), 


e a 23 
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It is now possible to find the reduction in power from 


Pa to Poe Applying the law of coSines to the diagram in 


Figure 5 resulte in 
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Normalizing, one obtains 
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where 

Е AR (C-2lı) 
is the per unit armature reactance loss. If rq. (C-23) 


1s solved Tor А then Ec. (11-15) can be evaluated 
oe a. 
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Ihe values ot Р thus obtained is compared with the 
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is incroased (or decreased) by the percentage difference 


between actual and desired output power (see statement 30 


in the design program). The program then returns to Eq. 


(C-18) or Eg. (C-19), and calculation continues in this 


fashion until an acceptable difference is reached. 


When a satisfactory P. is obtained the program continues 


with calculation of 
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The calculations described above are carried out for the 
initial value of Roe The program continues by incrementing 
Б, (by 0.1 in this program), returning to Eq. (C-11) and 
repeating all the steps through Eq. (€-208). Ro Continues 
to be incremented until it is found that tne volume for 
successive designs is increasing, which indicates that the 
minimum volume design is within the range of those designs 
that have been completed. The program then calls a sub- 
routine, described in Apvendix D, which evaluates the 
data for the designs that have been completed, and from 
this data determinss the minimum volume design. 

Within the desıen proceren, the following notation, 1s*used: 


РМА = desired Po 
POU e actual TP 


B(p,y) = orray of data representing the graphs 
in Figure 10 
OUTPT = array in which parameters of each design 
are stored 
B 
p max/const, 


XA = un-normalized armature reactance, which 
Decomes normalized wien suorea an Cult. 


т тт өз оңу r 
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ROS = R 
OS 


VOL = Volume 
МСТ = Veight 
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NIA = N Ia 


A print-out of the source dack for the design prosram 


3s included in tho following pages. 
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APPENDIX D 


Subroutine to Find the Minimum Volume Design 

After the design program computes three designs of 
successively increasing volume, it calls the subroutine 
MIN, The data stored in OUTPT et this time will include 


at least two designs beyond the minimum volume design. 
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The program computes dVol/dR, for the three values of 


R. that bracket tho minimum volume design. 
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A print-out of the source deck for this sub-routine is 


aociluded in the Lollwoings pages. 
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APPENDIX Ek 


Effoct of the Shield on Maximum Flux Density 

As described in Appendix B, the assumption was made, 
Mi determining characteristics of the fleld windings, that 
te effect of the shielding could be neglected. The 
Significance of this assumption is examined as follows: 

If Bo Zcoustoslseconpuucsuimas in Appendix P; Without 


meptecting shield effect, the result is 
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Assuming that ó is constant over a range of machine 
sizes, it is possible to compute the actual B p max/ Const. 
in terms of x and y. The solution technique is the seme 
as that described in Appendix B. One can then compare 
the value of B p max/ Const. obtained from the simplified 
equations with that value obtained from the method 


described above. 
The results of such a comparison are displayed in the 
following graphs. The data is for the cases of p = 3 to 6, 


which are of most interest. 
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1.2- Bo max With shield 
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Figure 14 


Bomax Without shield 
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Figure 15 


КОО Булак with shield 


В o max without shield 
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